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Self-consistency vanishes in the plateau regime of the bump-on-tail instability 
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Using the Vlasov-wave formalism, it is shown that self-consistency vanishes in the plateau regime 
of the bump-on-tail instability if the plateau is broad enough. This shows that, in contrast with 
the "turbulent trapping" Ansatz, a renormalization of the Landau growth rate or of the quasilinear 
diffusion coefficient is not necessarily related to the limit where the Landau growth time becomes 
large with respect to the time of spreading of the particle positions due to velocity diffusion. 
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I. INTRODUCTION 

The saturation of the bump-on-tail instability is a 
tough problem of kinetic plasma physics, which is still the 
subject of a controversy II, 2]. It was originally tackled in 
the frame of the Vlasov-Poisson formalism through the 
quasilinear approximation that neglects mode coupling 
[3, |4[ • However mode coupling was proved to be impor- 
tant during the saturation of the instability 5], which 
called for a new calculation scheme tailored to this phase. 
This is generically difficult because the particle motion is 
chaotic ; moreover wave-particle self-consistency makes 
the description of this chaotic regime especially difficult. 
This paper shows that, if the tail particle distribution 
displays a plateau, the waves whose phase velocity be- 
longs to the plateau feel a negligible mode-coupling in 
the limit where the plateau is broad enough. 



II. DESCRIPTION OF WAVE-PARTICLE 
SELF-CONSISTENCY 

The difficulty to describe the nonlinear regime of the 
Vlasov-Poisson system of equations, and the progress in 
the chaotic dynamics of Hamiltonian systems with a fi- 
nite number of degrees of freedom were an incentive to 
tackle the description of the saturation regime with the 
so-called self-consistent Hamiltonian that describes the 
one-dimensional self-consistent evolution, in a plasma 
with spatial periodicity L, of M Langmuir waves with N 
particles in the tail of the electron distribution function 
per length L [2[ ■ In the absence of tail particles, the Lang- 
muir waves are the collective motions of a plasma without 
resonant particles with a density ribuik, and with a plasma 
frequency u> p . Wave j has pulsation u>j related to the 
wavenumber kj (a multiple of 2n/ L) through the Bohm- 
Gross dispersion relation e(kj,u>j) = for the plasma 
dielectric function e. Each wave may be viewed as a har- 
monic oscillator. The interaction of these waves with the 
N tail particles is described by the self-consistent Hamil- 
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where (x r ,p r ) are the conjugate position and momentum 
of particle r (with mass normalized to unity), (Xj, Yj) are 
the conjugate generalized coordinate and momentum of 
the harmonic oscillator corresponding to Langmuir wave 
number j, /3j = [de/duj(kj, Wj)] _ > and 
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is the coupling coefficient, where m is the electron mass 
and ntaii is the density of the tail particles. This Hamil- 
tonian generates the evolution equations 

(3) 
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where K is the real part. Equation ((4]) makes clear the 
link between Zj = Xj + \Yj and the electric field of wave 

3- 

For the present paper, it is easier to describe the 
tail particles through a velocity distribution function 
f(x,p,t). This is possible through the so-called Vlasov- 
wave model that is obtained as a mean-field limit (limit 
N — > oo) of the dynamics defined by H sc @, 0], in anal- 
ogy to that performed by Spohn for his elegant and short 
derivation of the Vlasov equation j8|] ■ This yields the fol- 
lowing set of coupled equations 

M 

dtf + pd x f + n(j2Qe i(kiX ~" jt) ) d pf = 0. ( 6 ) 
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f(p, x, t)Apdx, 



(7) 



where A = [0, L] x R is the one-particle (x,p) space, and 
with the rescaled wave complex envelope variable 



Q = iepjZjJ' 



and coupling constant 



2fc; 
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Here J A f(p, x, t)dpdx = 1. This model keeps some sym- 
metry in the description of waves and particles, but gets 
rid of the granularity effects due to the discrete descrip- 
tion of particles, like spontaneous emission of waves by 
particles. 



III. DYNAMICS WHEN THE DISTRIBUTION 
IS A PLATEAU 

This paper considers the dynamics defined by Eqs (O 
[7]), while starting at time t = from (i) a spectrum of 
Langmuir waves where all nearby waves are in resonance 
overlap and (ii) a spatially uniform plateau for the par- 
ticle velocity distribution function over this overlap do- 
main with a height /o- We first start with a simplistic 
description of this dynamics, which will be useful to de- 
rive a more accurate one hereafter: since the plateau is 
spatially uniform, there is no source term for the waves in 
Eq. 0, which keep constant complex amplitudes; there- 
fore the particle dynamics is that defined by a prescribed 
spectrum of waves, which preserves the initial plateau. 
Clumps of particles may experience a strong turbulent 
trapping, but the distribution function is unaffected by 
this granular effect. 

In reality Kolmogorov-Arnold-Moser (KAM) tori 
bounding the chaotic domain defined by a prescribed 
spectrum of waves experience a sloshing motion due to 
the waves. This brings a small spatial modulation to 
the particle density which provides a source term for the 
Langmuir waves in Eq. IjTj). However the evolution of 
the wave spectrum is slow, which only brings a small 
change to the previous simplistic picture. This slow evo- 
lution suggests to introduce an adiabatic description of 
the true dynamics. 

To be specific, we consider the case where the wave 
phase velocities range over an interval [i>o,i>i], the near- 
est KAM tori have velocities v a and Vb and the particle 
distribution function plateau ranges over [uo,ui], with 
u o ^, v n ^ v o < v i ^S v h ^S u±. Therefore the plateau 
width Aw = Mi — wo is essentially equal to the chaotic 
domain width Vb — v a and to the wave spectrum width 
U\ — uq. Moreover, the wave phase numbers are all of 
comparable order of magnitude, typically k. Particles 
with velocity v experience the oscillations of a wave j 



at relative frequencies Qj 



— kjV, and for a parti- 



cle or a wave with velocity v the nearest KAM velocity 
defines f2 m i n = k min(t>b — v, v — v a ). The resonance over- 
lap in the wave spectrum is characterized by the small 



parameter k 1 / 2 5Q./C,J -C 1, where Sfl is the (Doppler) 
frequency detuning of two nearby waves and £* is the 
typical modulus of a wave amplitude. 

We now formalize this adiabatic description. Let (j (i) 
be the value of Q at time t in the self-consistent dynam- 
ics. Consider the non self-consistent dynamics D(io) de- 
fined by the self-consistent spectrum of Langmuir waves 
frozen at time to, as defined by Eq. ([6]) where the Q's are 
substituted with the Q(to)'s. The chaotic domain C(io) 
in single-particle (Boltzmann or p) phase space (x,p) de- 
fined by this frozen wave spectrum is bounded above 
and below in p by two KAM tori, respectively T a (to) 
and Tb(to)- The initial particle distribution function 
f(x,p,0) is assumed to be uniform on C(0) ; let /o be 
this uniform value. During the adiabatic evolution corre- 
sponding to the true self-consistent dynamics, f(x,p,t) 
stays uniform on C(i) and keeps the value /o- We are left 
with the calculation of the modulation of the width of a 
single water bag with height /o- The modulation of this 
width is given by that of KAM tori T a (t ) and T b (to), 
which may be computed by perturbation theory in the 
typical amplitude £* of the Langmuir waves. 

Since mode-mode coupling is a four-wave process, the 
first non vanishing contribution to Eq. Q is of order 
Q?. For waves with phase velocities further than Q m - ln /k 
from the edges of the plateau, a pessimistic estimate for 
(j(t) turns out to scale like (*Cl m i n /(AvSCl). With the 

(also pessimistic) estimate f2 m i n oc (* associated with 
the modulation of KAM tori, (* scales like ^ /(AvSQ), 
for a given spectrum with fixed discretization <5f2, and 
thus vanishes if the plateau is broad enough (Av — > oo). 
This justifies a posteriori our previous adiabatic approx- 
imation. As a result the plateau dynamics (further than 
^min/k from its boundaries) is almost the same as in a 
prescribed field of Langmuir waves. Therefore the chaotic 
motion of particles is almost unchanged due to the non- 
linear coupling of Langmuir waves. 



IV. CONCLUSION 

We have just shown that self-consistency vanishes in 
the plateau regime of the bump-on-tail instability if the 
plateau is broad enough. This means that the diffusion 
coefficient D(p) of particles with momentum p is that 
found for the dynamics of particles in a prescribed spec- 
trum of Langmuir waves. Let -Dql(p) be the quasilinear 
value of this coefficient. In the resonance overlap reg ime 
D/Dql may cover a large range of values @, M, [l0|- I n 
particular D ~ -Dql is obtained for random phases of the 
waves and strong resonance overlap 0, @, LLQ| • 

Let 7Landau be the maximum Landau growth rate 
of Langmuir waves for the instantaneous value of /. 
Let T sprca d = (/c 2 -D)~ 1//3 be the time after which the 
ballistic approximation fails for particles diffusing with 
the maximum instantaneous D due to a spectrum of 
Langmuir waves with typical wavenumber k. The 
plateau regime corresponds to 7LandauT spr cad = 0. Since 



D/Dql may cover a large range of values in this regime, 
7LandauT S prcad < 1 does not imply any renormalization 
of D/Dql, nor of 7/7Landau by wave-particle momentum 
conservation, in contrast with the "turbulent trapping" 
Ansatz of [ll| . The value of D/Dql in the plateau regime 



of the bump-on-tail instability depends on the kind of 
wave spectrum the beam-plasma system reaches during 
the saturation phase of the instability. A numerical sim- 
ulation using the Vlasov-wave model is now attempting 
to uncover this spectrum. 
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